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We discuss the prospects for setting limits on or discovering spin-1 Z' bosons using early LHC data 
at 7 TeV. Our results are based on the narrow width approximation in which the leptonic Drell-Yan 
Z' boson production cross-section only depends on the Z' boson mass together with two parameters 
c u and Cd- We carefully discuss the experimental cuts that should be applied and tabulate the 
theoretical next-to-next-to-leading order corrections which must be included. Using these results the 
approach then provides a safe, convenient and unbiased way of comparing experiment to theoretical 
models which avoids any built-in model dependent assumptions. We apply the method to three 
classes of perturbative Z' boson benchmark models: Eg models, left-right symmetric models and 
sequential standard models. We generalise each class of model in terms of mixing angles which 
continuously parametrize linear combinations of pairs of generators and lead to distinctive orbits in 
the c u — Cd plane. We also apply this method to the strongly coupled four-site benchmark model 
in which two Z' bosons are predicted. By comparing the experimental limits or discovery bands 
to the theoretical predictions on the c u -Cd plane, we show that the LHC at 7 TeV with integrated 
luminosity of 500 pb _1 will greatly improve on current Tevatron mass limits for the benchmark 
models. If a Z' is discovered our results show that measurement of the mass and cross-section will 
provide a powerful discriminator between the benchmark models using this approach. 

I. INTRODUCTION 

The end of the first decade of the millenium is an exciting time in particle physics, with the CERN LHC enjoying an 
extended run at 7 TeV, and the Fcrmilab Tevatron collecting unprecedented levels of integrated luminosity, eventually 
up to perhaps 10 fb _1 , in the race to discover the first signs of new physics Beyond the Standard Model (BSM). Since 
spin-1 Z 1 bosons are predicted by dozens of such models, and are very easy to discover in the leptonic Drell-Yan 
mode, this makes them good candidates for an early discovery at the LHC. For a review see Refs. [l|-[3| and references 
therein. Furthermore high mass Z 1 bosons are more likely to be discovered at the LHC than the Tevatron since 
energy is more important than luminosity for the discovery of high mass states. This makes the study of Z' bosons 
both timely and promising and has led to widespread recent interest in this subject (see for example |5H2ll l2l|). 

Since one of the purposes of this paper is to facilitate the connection between experiment and theory, it is worth 
being clear at the outset precisely what we shall mean by a Z' boson. To an experimentalist a Z' is a resonance 
"bump" more massive than the Z of the Standard Model (SM) which can be observed in Drell-Yan production 
followed by its decay into lepton-antilepton pairs. To a phcnomenologist a Z 1 boson is a new massive electrically 
neutral, colourless boson (equal to its own antiparticle) which couples to SM matter. To a theorist it is useful to 
classify the Z' according to its spin, even though actually measuring its spin will require high statistics. For example 
a spin-0 particle could correspond to a sneutrino in R-parity violating supersymmetric (SUSY) models. A spin-2 
resonance could be identified as a Kaluza-Klein (KK) excited graviton in Randall-Sundrum models. However a spin-1 
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Z' is by far the most common possibility usually considered, and this is what we shall mean by a Z 1 boson in this 
paper. 

In this paper, then, we shall discuss electrically neutral colourless spin-1 Z 1 bosons, which are produced by the 
Drcll-Yan mechanism and decay into lepton-antilepton pairs, yielding a resonance bump more massive than the Z . 
We shall be particularly interested in the prospects for discovering or setting limits on such Z' bosons using early 
LHC data. By early LHC data we mean the present 2010/11 run at the LHC at 7 TcV, which is anticipated to yield 
an integrated luminosity approximately of 1 fb _1 . Since the present LHC schedule involves a shut-down during 2012, 
followed by a restart in 2013, the early LHC data will provide the best information possible about Z 1 bosons over 
the next three years, so in this paper we shall focus exclusively on what can be achieved using these data, comparing 
the results with current Tevatron limits. In order to enable contact to be made between early LHC experimental 
data and theoretical models, we advocate the narrow width approximation, in which the leptonic Drcll-Yan Z' boson 
production cross-section only depends on the Z' boson mass together with two parameters c u and Cd Properly 
defined experimental information on the Z' boson cross-section may then be recast as limit or discovery contours in 
the c u — Cd plane, with a unique contour for each value of Z' boson mass. In order to illustrate how this formalism 
enables contact to be made with theoretical models we study three classes of Z' boson benchmark models: E§ models, 
left-right (LR) symmetric models and sequential standard models (SSM). We also apply this method to the strongly 
coupled four-site benchmark model in which two Z' bosons are predicted fl3j . Each benchmark model may be 
expressed in the c u — Cd plane which enables contact to be made with the experimental limit or discovery contours. 
Working to next-to-next-leading order (NNLO) wc show that the LHC at 7 TeV with as little data as 500 pb _1 
can either greatly improve on current Tevatron mass limits, or discover a Z', with a measurement of the mass and 
cross-section providing powerful resolving power between the benchmark models using this approach. Wc also briefly 
discuss the impact of the Z' boson width on search strategies. Although the width Tz' into standard model particles 
may in principle be predicted as a function of Mz> , in practice there may be considerable uncertainty concerning L ' z 1 
due for example to the possible decay into other non-standard particles, including supersymmetric (SUSY) partners 
and exotic states, for example. 

At the outset wc would like to highlight some of the strengths and limitations of our approach to the benchmark 
models. One of the strengths of our approach is that the considered benchmark models encompass two quite different 
types of Z' models: pcrturbativc gauge models and strongly coupled models, where the perturbative models generally 
involve relatively narrow widths (which however can get larger if SUSY and exotics are included in the decays in 
addition to SM particles), while the strongly coupled models involve multiple Z' bosons with rather broad widths. 
The perturbative benchmark gauge models are defined in terms of continuous mixing angles, in analogy to the Eq 
class of models expressed through the linear combinations of x an d ip generators. This approach is generalised to the 
case of LR models involving linear combinations of the generators R and B — L 1 , and similarly we define a new class 
of SM-likc Z' models involving linear combinations of L and Q generators, with the precise linear combination in each 
case parameterised by a separate angle. The strength of this approach is that it enables a finite number of classes, 
each containing an infinite number of benchmark models, to be defined, rather than just a finite number of models. 
For each class of models, the respective angle serves to parametrize the specific orbit which represents that class in the 
c u -Cd plane. The different orbits turn out to be non-overlapping for the abovementioned three classes of models which 
in the following we label as: Eq(9), GLR((f>) and GSM (a) respectively. A limitation of the approach is that it ignores 
the effect of the SUSY and exotics on the width L 'z' ■ In addition it also ignores the effects of Z — Z 1 mixing since 
this is model dependent. However any such mixing must be small due to the constraints from electroweak precision 
measurements, and we refer to such constraints on the mixing angle where possible. As regards the strongly coupled 
models, in principle Z' bosons could emerge from techni-rho bound states in Walking Technicolour models, or a scries 
of strongly interacting resonances such KK excitations of the Z boson. A limitation of our approach here is that we 
only consider the four-site Higgless model which contains just two excitations of the Z (and W) bosons, as a simple 
approximation to both the Walking technicolour models and the KK excitations of the Z. Neverthless the four-site 
model is representative of the physics of a typical strongly interacting Z' model and by representing it for the first 
time in the c u -Cd plane, it is clearly seen that the associated Z' bosons may easily be distinguished from those of the 
perturbative gauge models. 

The layout of the rest of the paper is as follows. In section [II] wc describe our model independent approach based on 
the narrow width approximation and the variables Mz>, c u and Cd- Higher order corrections to the cross-section are 
tabulated in the narrow width approximation and new K factors are defined which enable the c u and Cd approach to 



1 We remark that in |19| the authors focussed on models obtained by taking a linear combination of Y and B — L which is related to the 
generalized LR models that we also consider, however the choice of generator basis and paramctrisation is different and they did not 
use the mixing angle paramctrisation that we propose here. 
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be reliably extended to NNLO. In this section we consider finite width effects and discuss the choice of invariant mass 
window around Mz> which matches the narrow width approximation, showing that in the case of the LHC this choice 
is only weakly constrained. We also comment on the effects of interference and show that they may become important 
for invariant mass cuts close to 100 GcV but are negligible for the a suitable invariant mass cut around Mz>, which 
we therefore advocate. In section Hill we define our benchmark models based on Eq, LR and SSM, generalised using 
variables which continuously parametrize mixing of the respective f (1) generators. 

For these perturbative models we specify the gauge coupling and calculate the vector and axial couplings in terms of 
the mixing angle variable. We tabulate the results for some special choices of the mixing angle variable which reproduce 
the models commonly quoted and analysed in the literature. For the four-site model we discuss the parameter space 
describing the two Z' bosons allowed by electroweak precision measurements and unitarity. In section IIVI we apply 
the model independent approach to the benchmark models discussing first the current Tevatron limits using the latest 
DO results with collected luminosity L = 5.4/6 -1 , then the expected LHC potential based on the projected CMS 
limits on the Z' boson cross section normalized to the SM Z-boson one for an integrated luminosity L = 500pb~ 1 . 
We use the DO results rather than CDF since they are more closely related to the narrow width approximation that 
we advocate, and we use CMS rather than ATLAS since the projected limits are publicly available. In both cases we 
express the experimental cross-section limits in the c u -Cd plane and compare these limits to the benchmark models 
also displayed in the c u -Cd plane. In the case of CMS we also show the discovery limits. We tabulate some of the 
results for some special choices of the mixing angle variable, including the width, the leptonic branching ratio, c u and 
Cd values, our derived current direct limit on the Z' mass based on DO results, as well as other indirect and mixing 
limits where available. Finally in section IVl we discuss the impact of Tz> on search strategies. Section |VT1 summarises 
and concludes the paper. 

II. MODEL INDEPENDENT APPROACH 
A. Couplings 

At collider energies, the gauge group of a typical model predicting a single extra Z' boson is: 

SU(3) C x SU(2) L x U(1) Y x U'(l) (II.l) 

where the Standard Model is augmented by an additional U'(l) gauge group. The U'(\) gauge group is broken near 
the TeV scale giving rise to a massive Z' gauge boson with couplings to a SM fcrmion / given by: 

The values of g v , g A depend on the particular choice of U'(l) and on the particular fermion /. We assume universality 
amongst the three families, g v = g v = g v , and gy = gy — gy, as well as g v = g v = g v , and similarly for the g A 
couplings, which means that there are eight model dependent couplings to SM fermions g v A , with / = u,d,e,v. 
These eight couplings are not all independent since they are related to the left (L) and right (R) couplings as follows: 

£nc = ^Z'J-ftfv - .^7 5 )/ = g'Z'J^(e{P L + e f R P R )f. 

where Pl.r = (1 T 7s)/2, and g v A = e£ ± e^,. The couplings are not all independent since the left-handed fermions 
are in doublets with the same charges = t d L and e e L = e u L . Excluding the right-handed neutrinos (which we assume 
to be heavier than the Z') there are really five independent couplings ef, =1/ , ^t =d , efj- However we prefer to work 

with the eight vector and axial couplings g{/ A - In addition, the strength of the gauge coupling g' is model dependent. 
Throughout, we follow the conventions of [l|. 

A slightly more complicated setup is needed to describe the four-site model which, in this paper, has been chosen 
to represent Higgslcss multiple Z'-boson theories. The corresponding framework will be given in Scc lIII B~T1 

Throughout the paper we shall ignore the couplings of the Z' to beyond SM particles such as right-handed neutrinos, 
SUSY partners and any exotic fermions in the theory, which all together may increase the width of the Z' by up to 
about a factor of five [22J and hence lower the branching ratio into leptons by the same factor. 
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B. Z 1 production and decay in the narrow width approximation 



The Z 1 contribution to the Drell-Yan production cross-section of fermion-antifcrmion pairs in a symmetric mass 
window around the Z' mass (|M — M%t\ < A) may be written as: 



,(M Z ,+A) 2 ■ 

r n = ^7T2(PP^ Z '^ ff X )dM 2 . {11.2) 

11 J(M Z ,-Ay dM 



In the narrow width approximation (NWA) , it becomes 



°fj « f \ E (§jf) "> ^) J >< MZ' -> //) (IL3) 

V q=u,d \ Z' / I 

where the parton luminosities are written as {^jjfr-^ arid cr^g — > Z') is the peak cross-section given by: 

a{ q q^Z') = y\glf + {g\f]. (II.4) 
The branching ratio of the Z' boson into fermion-antifermion pairs is 

Br(ff) ee Br(Z' -+ //) = T{Z ' ^ //} (II.5) 

where Tz> is the total Z' width and the partial widths into a particular fermion-antifcrmion pair of N c colours is given 
by 

T(Z' H- //) = N c ^-M z ,[{gl? + (gt A )% (11.6) 

Assuming only SM fermions in the final state and neglecting in first approximation their mass, one finds the total 
width: 



Tz ' = W, Mz ' 



H9v 2 +9a)+H9v +9A) + H9v 2 +9A 2 )+3(g"v 2 +9l 2 ) ■ (H.7) 



Specializing to the charged lepton pairproduction cross-section relevant for the first runs at the LHC, Eq III. 31 may 
be written at the leading order (LO) as [J]: 

(Ti+ e - = ^ [c u w u (s, M%) + c d w d (s, M%)\ (II.8) 

where the coefficients c u and c d are given by: 

c u = 9 ^-(9 u v 2 + 9 U A 2 )Br(£+r), c d = 9 -^(gf + g d A 2 )Br(i+ T\ (II.9) 

and w u (s, M§,) and w d (s, Mf,) are related to the parton luminosities ( ^ and {^jff-^j and therefore only depend 

on the collider energy and the Z' mass. All the model dependence of the cross-section is therefore contained in the 
two coefficients, c u and c d . These parameters can be calculated from gy,g A and g' , assuming only SM decays of the 
Z' boson. The corresponding values for all models which predict a single Z' boson purely decaying into SM fermions 
are given in Table [TTJ 

A slight complication arises in Higgsless theories, which in the present paper are represented by the four-site model. 
Here in fact the two neutral extra gauge bosons, Z1.2, decay preferebly into di-boson intermediate states. Their total 
width gets therefore two contributions: 

rz 4 =r£{ + i% v (i = i,2) (11.10) 

where the two terms on the right-hand side represent the fermionic and bosonic decay, respectively. More in detail, 
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9(<7,V 2 + 9Xa 2 ) + H9tv 2 + <£a) + 3(5fy 2 + qIa) + Kd^ + 9^)} (H.ll) 
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r "' w 'i{i) 2 W {1 - z,){1+z2) (IU2) 

with i=l,2, where in this case we have included the coupling in the definition of g( 2V and g{ 2 A- ^ n the above 
formulas Mzi and Mzi are the masses of the two extra gauge bosons, Z\,2, while z is their ratio, i.e. z = Mzi/Mz2- 
The direct consequence of this peculiarity is that the Z\ y 2 leptonic branching ratio acquires a not trivial dependence 
on the Z12 boson mass which reflects in an intrinsic mass dependence of the c u and Cd coefficients. In addition, 
there is an external source of variation with mass. As all vector and axial couplings in the four-site model can be 
expressed in terms of the three independent free parameters (3^, Mz2, z), c u and c c i are completely specified by 
these quantities as well: c u _d = Cti.cKfffyi Mz2, z). This means that at fixed masses, Mz2 and z = Mzi/Mz2, these 
coefficients get constrained by the EWPT bounds acting on g 2V . As these limits vary with mass (see Fig. [6]), c u and 
Cd acquire this extra Mz2 dependence. The net result opens up a parameter space in the Cd — c u plane which will be 
dispalyed at due time. 

As emphasized in Q , the Cd — c u plane parametrization is a model- independent way to create a direct correspondance 
between the experimental bounds on pp(p) — > Z' — > £ + £~ cross sections and the parameters of the Lagrangian. An 
experimental limit on a(pp(j)) — > Z' — > £ + £~) for a given Z' mass gives in fact a linear relation between c u and Cd, 

c u =a-bc d (H.14) 

where a, b can be regarded as known numbers given by: 



7T W u 



cxp 

b=^. (11.15) 



A8sa e+e _ t w d 



where cr^+^- represents the 95% C.L. upper bound on the experimental Drell-Yan cross section which can be derived 
from observed data. 

In practice, it is more convenient to use a log-log scale resulting in the limits appearing as contours for a fixed Z' 
mass in the Cd ~ c u plane. We use this representation in the next subsections. 



C. Higher-order corrections 

At higher-orders, the expression for Z' production given by Eq. (|II.8|) strictly speaking is no longer valid. However, 
as it was shown in Ref. [3] , the additional terms which are not proportional to c u and Cd in Eq. (|II.8|) can be neglected 
at NNLO. Therefore, Eq. (|II.8[) gives a quite accurate description of the approach we are discussing here even at 
NNLO. 

In the following, we take into account QCD NNLO effects as implemented in the WZPROD program p3l425| as a 
correction to the total Z' production cross section in the NWA. 2 We have adopted this package for simulating the Z' 
production, and have linked it to an updated set of Parton Density Functions (PDF). This set includes in particular 



the most recent versions of CTEQ6.6 [29j, |30( and MSTW08 [31( PDF's, which we use in our analysis. We can provide 
the complete code upon request. 

The QCD NNLO Z' production cross sections obtained using CTEQ6.6 and MSTW08 PDFs are in a good agree- 
ment. Their difference is in fact at the 2-3% level over a wide Z' mass spectrum as shown in Fig. [TJ where we plot 
the total pp(p) — > Z' cross section at the Tevatron and LHC@7TeV versus the Z' mass. Here, we have taken as 
factorization scale the value Q = Mz>- The further detailed analysis of the cross section variation with the scale is 
outside of the scope of the current paper. 

It is also convenient to define customary NLO and NNLO 7^-factors which can be useful for experimentalists in 
establishing Z' exclusion limits: 

cr(pp(p) -> Z')j 

Ki = (11.16) 

cT{pp(p) -> Z')o, 



2 We would like to note, that study of NLO effects for kinematical distributions involving leptons from Z' decay I26H28H is beyond the 
scope of this paper. 
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FIG. 1: cr(pp(p) -¥ Z') NN lo for Standard Model-like Z' production at the Tevatron (left panel) and the LHC@7TeV (right 
panel) for CTEQ6.6 and MSTW08 PDF's 
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200 400 600 800 1000 1200 1400 500 1000 1500 2000 2500 3000 
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FIG. 2: NNLO and NLO K -factors defined by Eq. (|IL16|l for Standard Model- like Z' production at the Tevatron (left panel) 
and the LHC@7TeV (right panel). 



where the index i = 1,2 corresponds to NLO and NNLO K- factors, respectively. As an example, in Fig. [2] we present 
the values of these JQ-factors for Standard Model-like Z' production at the Tevatron (left panel) and the LHC@7TcV 
(right panel) for CTEQ6.6 and MSTW08 PDF's. 

Oppositely to what happens for the aforementioned exact NNLO Z' production cross section, where the agreement 
between CTEQ6.6 and MSTW08 PDF predictions is optimal, in this case there is a noticablc difference in the behavior 
of Knlo and Knnlo factors as a function of the Z' mass when convoluting the Z 1 production cross section with 
CTEQ6.6 or MSTW08 PDF's. This difference is related to the way of fitting the LO PDF's of CTEQ and MSTW 
collaborations (see e.g. [3l|,[32|). Furthemore, both Knlo^n'nlo an< ^ ^nlo^nnlo iac t° rs display a strong dependence 
on the Z' mass. As an example, -Kjvwlo 6 var i es between 10-40% at the Tevatron and 10-30% at the LHC@7TeV for 
potentially accessible Z' masses. 

Applying a universal X-factor can be highly misleading. As shown above, the Knlo,nnlo-£&c^ot has indeed a 
two- fold source of dependence: PDF set and energy scale (i.e. Mz>)- A uniform setup must be fixed when comparing 
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experimental limits on different models. 

Since Eq. (|II.8[) gives an accurate description even at NNLO and noting that QCD NNLO corrections are 
universal for up- and down-quarks, one can effectively apply the same K jv n xo-factor derived for SM-likc Z' to 
generic Z' models without loosing of generality Owing to the remarkable Z' mass dependence of the K n n LO-fa c t° r , 
we first convolute the LO Z' production cross section with the respective LO PDF's and then we multiply it by 
K NNLO (M Z '). 

For convenience and clarity, we provide in Tables Hill and IIVI shown in Appendix lAl the values of K jy n Lo _ factors 
and cross sections for the SM-like Z'-boson production process at the Tevatron and the LHC@7TeV: pp(p) — > Z' + X. 
The first table contains the results obtained with MSTW08 PDF, the latter with CTEQ6.6 PDF. The quoted numbers 
correspond to the curves visualised in Figs lll2l 

In narrow width approximation, the two-fermion cross section is the product of the production cross section and the 
respective branching ratio. When considering the complete Z'-boson production and decay in the Drcll-Yan channel, 
one has to keep in mind that QCD NNLO corrections also affect the Z' branching ratio even for purely leptonic 
decays, Br(Z' —> £ + £~), since the Z' total decay width will be corrected at NNLO. This reflects into an higher order 
correction to the c u and c d coefficients, through Br(Z' — > £ + £~) which explicitly enters the expression for c u and c d 
given in Eq. pi.9[) . The NNLO Drell-Yan cross section can be thus written as: 

NNLO n r NNLO i A r 2 \NNLO , NNLO ( A j2 \NNLO~l 

e+i ~ = 48s ^ u w u{s,M z ,) +c d w d (s,M z ,) \ 

= K% N F LO K% NLO -^- [c u w u (s, Ml) + c d w d (s, M%,)] = K™ F LO K™ lo ^°- (11.17) 

The leading NLO QCD correction to the total Z' width is known to be a s /Tr [3^, [H|. This gives an enhancement 
of the order of 2-3% to the Z' width for M z > in the range 500-2000 GeV. The Br(Z' -> £+£~) will thus decrease 
accordingly by (2 — 3%) x Br(Z' — > hadrons). The net result corresponds to a 1 — 2% deplection of the leptonic 
branching ratio within the SM-likc Z' model. An effect of the same order is expected for the other classes of Z' models 
under consideration. In the current study, we neglect this effect and use the following formula for establishing limits 
on Z' models: 

NNLO ^ t,-PDF LO m i o\ 

V>t+e- — K NNLO a t+t-- (11.18) 



D. Finite width effects 



So far we have have discussed the Z' boson production using narrow width approximation. However, the exper- 
imental search for an extra Z' boson and the descrimination of the SM backgrounds could strongly depend on the 
realistic Z' width. Moreover the theoretical prediction of the Z 1 production cross section also depends on its width 
as we discuss below. 

We start this discussion with Fig. [3] where we present the di-lepton invariant mass distribution for the Z' boson 
production within various models at the Tevatron (left panel) and the LHC@7TeV (right panel). We consider three 
representative models: the SM-like Z 1 model (black line), the N-type Eq model defined in Table Q] (red line), and the 
weakly coupled SM-like Z' model where the Z' boson gauge coupling to SM fermions is reduced by a factor 10 (blue 
line). From top to bottom, the last two distributions are normalized to the integral under the first one. We first 
consider the SM-like Z 1 model distribution at the Tevatron. It is important to stress that the total cross section of 
pp —> Z' — > £ + £~ process integrated over the entire Mg+£- range is actually is almost as twice as large as the SM-like 
Z' in the narrow width approximation. The main reason for this effect is the specific shape of the M e + e - distribution 
in the region of small Mg+g- far away from M Z i. This region is exhibited by a non-negligible tail due to the steeply 
rising PDF in the region of low M t + g- even though the Z' boson isextremely far off mass-shell in this region. The 
integral over this region can even double the cross section evaluated in the NWA in the case of Z' production at the 
Tevatron. 

This effect, which is related to the off-shellness of the extra gauge boson, varies according to the total Z' width. 
In the Z' N model, it brings an additional 20% contribution to the narrow width approximation cross section at the 
Tevatron. In the weakly coupled SM-likc Z' model, the far off-shellness effects are effectively negligible (below 1%). 

We can see that in general experimental limits would and should strongly depend on the particular Z' model 
predicting a specific Z' width. On the other hand, if one requires a di-lepton mass window cut around the Z' mass, 
one can establish a quasi model-independent experimental upper limit on a(pp — > Z' x Br(Z' — > £ + £~) versus M z > 
and apply this limit to constraint different classes of models. 

In Fig. U we present the effect of a symmetric mass window cut around M z > for the SM-likc Z' model and two other 
representative models (see Table fl} at the Tevatron and the LHC@7TcV. We fix the Z' mass to be M Z i = 1 TeV. 
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FIG. 4: Relative difference between the full cross section for pp(p) — > Z' — > £ + £~ evaluated taking into account the finite Z' 
width (a) and the cross section computed in narrow width approximation (a NWA ). The relative difference is presented as 
a function of the AM/M Z ' symmetric mass window cut \M e + £ - — M Z '\ < AM applied to the full cross section (<r). Three 
different representative models are considered for M z , = 1 TeV (see Table [TJ. 



We plot the relative difference between the full cross section for the process pp(p) — > Z' — > £ + £~ evaluated taking 
into account the finite Z' width (a) and the cross section computed in narrow width approximation (a NWA ). The 
relative difference is presented as a function if the AM/Mz' symmetric mass window cut (\Mg+£- — Mz*\ < AM) 
applied to the full cross section a. One can see that for the SM-like Z' model at the Tevatron, a AM/M cut in the 
9-25% range brings the agreement between a and a NWA down to 5% level while AM/M ~ 15% exactly matches a 
and cr NWA . At the LHC, the corresponding range of the AM/M cut is 15-80%, and a and a NWA are matched for 
AM/M ~ 45%. The if: model Z' has a narrower width, making the choice of the cut more insensitive, while the Q 
model Z' width is broader leading to a more sensitive choice of the mass window cut to reproduce the narrow width 
approximation. Note that all lines cross the abscissa at about the same value of AM/Mz>, meaning that there will 
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be an optimal mass window cut consistent with all models. The choice of the mass window cut to gain agreement 
with the narrow width approximation also depends on Mz> ■ This dependence is defined by proton parton densities 
and is therefore model- independent. The net effect is again to make all the lines cross the abscissa at about the same 
value of AM/Mz 1 , where this point depends on Mz> ■ Therefore for every given mass one can work out a quasi model- 
independent mass window cut where the full cross section matches the narrow width approximation. The additional 
advantage of this choice is that in the selected mass window around the Z' mass the model-dependent interference 
effect between Z' signal and SM background is highly suppressed. 

The experimental limits would be quasi model-independent if one would apply this cut on the Mi+g- around the 
M' z : it brings in agreement the cross section calculated in the narrow width approximation and in the finite width 
approximation as well as removes model-dependent shape of the Mg+g- distributions in the region of low Mg+g- 
cspecially for the case of large Z' width effects as, for example, take place for SM-like Z'. Moreover, the cut on M^+^~ 
around the M' z plays an important role in reducing an effect of Z' interference with Z/j down to a few% level, which 
again, allows to establish and use experimental limits in model-independent way. 

For example, in case of SM-like Z' production at the Tevatron, the relative interference, which is defined as 
Ri = [cr(pp -> Z'/Z/i -> t+t~) - <r(pp -> Z' -> £+£-) - a(pp -> Z/ 7 -> £+£-)}/a(pp -> Z' -> £+T) is as 
large as about —19 (meaning —1900% of interference^)) for M e + e - > 100 GeV cut but it drops down to —6% for 
\Mg+£- — Mz'\ < 0.15Mz> cut which matches NWA and finite width cross sections. The effect of the mass window 
cuts is also quite large for the case of SM-like Z' production at the LHC, where interference is about —300% for 
M i+i - > 100 GeV cut and only about -2% for \M i+g - - M z > \ < 0.15M Z ' cut. 

We can see, that there is a strong motivation to use an invariant mass window cut for conducting a model- 
independent analysis. The size of this cut, if one aims to match the NWA and finite width cross sections, is collider 
dependent: it is about 15% of M z > for the Tevatron and about 40% of M z < for the LHC7TeV. 

In this paper we are using results of experimental analysis which are based on M e +£- mass window cut similar 
to what we are advocating. This would allow us to use precise NNLO model predictions and perform a respective 
model-independent interpretation of the experimental limits. 

III. BENCHMARK MODELS 

In this section, we extend and classify the benchmark Z' models present in the literature. We divide such classes 
into two main types: perturbative and strongly coupled gauge theories. 

A. Perturbative gauge theories 

1. E 6 Models 

In these models one envisages that at the GUT scale the gauge group is Eq, The gauge group Eq is broken at the 
GUT scale to 5*0(10) and a U(l)^ gauge group, 

E 6 -> 50(10) x U(l)^. (III.l) 

The 5O(10) is further broken at the GUT scale to SU(5) and a U(l) x gauge group, 

5O(10) -> 577(5) x U(l) x . (111.2) 

Finally the SU (5) is broken at the GUT scale to the Standard Model (SM) gauge group, 

SU(5) -> SU(3) C x SU(2) L x U(1) Y - (111.3) 

All these breakings may occur at roughly the GUT scale. The question which concerns us here is what happens to 
the two Abelian gauge groups U(l)^, and U(l) x with corresponding generators T^, and T x . Do they both get broken 
also at the GUT scale, or may one or other of them survive down to the TeV scale? In general it is possible for some 
linear combination of the two to survive down to the TeV scale, 

u(iy = cose u(i) x + sine u(i)^, (111.4) 

where —tt/2 < 6 < tt/2. More correctly the surviving Eq generator Qe 6 should be written as, 

Q Ee = cos 0T X + sin 6 T 4 , . (III. 5) 



10 



Some popular examples of such U(l)' are shown in Table |U 

The resulting heavy Z' couples as g'QE 6 Z'. Note that in E§ models it is reasonable to assume that the Z 1 gauge 
coupling g' is equal to the GUT normalized U(l)y gauge coupling of the SM, g\{Mz) = (e/cw)y/5/3 ~ 0.462 where 
e = 0.3122(2) and c w = y/l - s 2 w where the 175 value is s\ = 0.2312. Thus we take g' « 0.46. GUT normalization 
also implies that the charges of the fermions in the 5O(10) 16 representation for the tp case are all equal to l/y/24, 
while for the x case the T x charges of the SU(5) representations (10, 5, 1) arc (— 1/V40, 3/V40, — 5/v40). Recalling 
that g{r a = e{ ± £p, and 10 — > Q,u c ,e c and 5 — > L,d c , and that u c ,d c ,e c have the opposite charges to u R ,d R ,e R , 

this results in the values of the g(r A charges for the U(l)^ and U(l) x cases as shown in Table HI The general charges 
as a function of 8 are then simply given as, 

g f VA (8) = cosd g f VA ( X ) + sin0 g f v>A (i>), (III.6) 
where the numerical charges for the popular models quoted in the literature are listed in Table HI 



2. Generalised Left-Right Symmetric Models (GLRs) 

These models are motivated by the left-right (LR) extensions of the SM gauge group with the symmetry breaking, 

SU{2) L x SU{2) R x U(1) B -l -> SU{2) L x U{l) Y (III.7) 

which, from the point of view of Z' models essentially involves the symmetry breaking, 

U(1) R x U{1) B -l ->■ U(1) Y (III.8) 

where U(1) R involves the generator T 3R corresponding to the third component of SU(2) R , while U(1)b-l involves 
the generator T B -l = (B — L)/2. The hypercharge generator is then just given by Y = T 3R + T B _ L . Assuming 
a left-right symmetry, the gauge couplings of SU(2)l,r are then equal, gi, = g R and the resulting heavy Z' LR then 
couples as giQhRZ' where 




(III.9) 



with a = \/ cot 2 0w — 1 ~ 1-53 and g\ « 0.462 as before. 

The left-right symmetric models therefore motivate a U(1)lr which is a particular linear combination of U(1) R 
and U(1)b—l with a specific gauge coupling. From this perspective the special cases where the Z' corresponds to 
a pure U(1) R or a pure U{\)b-l ar e not well motivated. Nevertheless these types of Z' have been well studied in 
the literature and so it is useful to propose a generalization of the LR models which includes these special cases. To 
this end we propose a generalized left-right (GLR) symmetric model in which the Z corresponds to a general linear 
combination of the generators of U(1) R and U(1)b-l, 

Qglr = cos cj) T 3R + sin T B _ L , (III. 10) 

where —ir/2 < <f> < ir/2. The gauge coupling g' is fixed so that for a particular value of <j> the Z' of the GLR may be 
identified with the Z' of the LR symmetric model above. To be precise, we identify, for a particular value of <j): 

,9iQlh = g' Qglr (HI.ll) 

which implies tan0 = — 1/a 2 which corresponds to <j> = — 0.1287T for a w 1.53 and we find g' = 0.595. Keeping 
g' = 0.595 fixed, we are then free to vary </> over its range where <j) = — 0.1287T gives the LR model, but other values 
of <j) define new models. 

Clearly = gives a U(1) R model while <f> = tt/2 gives a U(1)b-l model. In the GLR model the value of <j> = 7r/4 
also defines a Z' which couples to hypercharge Y = T^r+Tb-l ( n °t to be confused with the sequential SM Z 1 which 
couples like the Z) . The couplings of the Z 1 for the special cases of the GLR models are give in Table HI The general 
charges as a function of <fr are then simply given as, 

g^ A {<P) = co$<t>g^ A {R)+sm4>g^ A {B-L), (111.12) 
where the numerical charges for particular models are shown in Table HI 
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TABLE I: Benchmark model parameters and couplings. The angles 9, <f>, a are defined in the text. 

3. Generalised Sequential Models (GSMs) 

No study is complete without including the sequential standard model (SSM) Z' SSM which is defined to have 
identical couplings as for the usual Z, namely given by —QzZ' SSM and Qz = T?,l — s^Q where = 0.2312 and 
ui{Mz) = 5l/(47 r ) ~ 0.0338 imply that « 0.74. Similar to the GLR models, it is useful to define a generalised 
version of the SSM called GSM where the heavy gauge boson Z' GSM then couples as g' QgsmZ' gsm where Qgsm 
corresponds to a general linear combination of the generators of U(1)t 3L and C/(1)q, 

Qgsm = cos a T aL + sin a Q, (111.13) 

and where — ir/2 < a < ir/2. The gauge coupling g' is fixed so that for a particular value of a the Z' GSM of the GSM 
may be identified with the Z' SSM of the SSM above. To be precise, we identify, for a particular value of a: 

—Qz = g'QasM- (HI.14) 

c w 

This implies that the GSM reduces to the SSM case for g' = —^Jl + s^, rj 0.76 and tana = —0.23 which corresponds 
to a = — 0.0727T. Keeping g' = 0.76 fixed, we are then free to vary a over its range where </> = — 0.0727T gives the usual 
SSM, but other values of a define new models. Clearly a = gives a U(1)t 3L model while a = ir/2 gives a U{1)q 
model. 

The couplings of the Z 1 for the special cases of the GLR models are give in Table HI The general charges as a 
function of a are then simply given as, 

9v^( a ) = cosa 9v,a( l ) + sina 9v,a(Q)> (111.15) 
where the numerical charges for particular models are shown Table |TJ 

B. Strongly coupled gauge theories 

Strongly interacting gauge theories provide an alternative mechanism for the electroweak symmetry breaking 
(EWSB). The EWSB is not driven by a light Higgs boson anymore, but it happens in a dynamical way. Such 
theories date back to decades. However, even if they predict the existance of new gauge bosons in order to delay at 
high energy the perturbative unitarity violation in vector boson scattering amplitudes, they are not considered when 
performing searches of Z' bosons in the dilepton Drcll-Yan channel. The reason is that historically the predicted 
new resonances must be fermiophobic in order to evade EWPT constraints. However, in recent years, new models 
have been proposed that are able to satisfy the EWPT bounds without imposing such a strong condition. Both the 
Minimal Walking Technicolour [HI, HH and the four site Higgsless model fl3l [37l |3&| predict extra Z 1 bosons with 
sizeable couplings to SM matter. Hence, they could be tested in the favoured Drcll-Yan channel at the Tevatron and 
during the early stage of the LHC. 
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FIG. 5: Left: Unitarity bound as a function of the energy scale for different z = Mzi/Mz2 values. The perturbative region is 
on the left of the contour plot. Right: 95% C.L. bounds on the vector couplings of the Zi^-bosons to SM electrons from ei 
(blue) and £3 (green). We consider the representative case: Mzi = 1 TeV, Mz2 — 1.3 TeV. 



1. The Four Site Higgsless Model 



Higgslcss models emerge naturally from local gauge theories in five dimensions. Their major outcome is delaying 
the unitarity violation of vector-boson scattering (VBS) amplitudes to higher energies, compared to the SM without 
a light Higgs, by the exchange of Kaluza-Klein excitations [33 ]. Their common drawback is to reconcile unitarity with 
the Electro Weak Precision Test (EWPT) bounds. Within this framework, and in the attempt to solve this dichotomy, 
many models have been proposed [40l - [48| . 

In this paper, we consider the four site Higgsless model [49| as rep resentative of strongly coupled theories. This 
model belongs to the class of the so called deconstructed theories [50l - l58j which come out from the discretization of 
the fifth dimension on a lattice, and are described by chiral lagrangians with a number of gauge-group replicas equal 
to the number of lattice sites. The simplest version of this class of models is related to the old BESS model (5(J H(|, a 
lattice with only three sites and SU(2)l x SU(2) x U(1)y gauge symmetry (for it, sometimes called three-site Higgsless 
model). In order to reconcile unitarity and EWPT-bounds, this minimal version predicts indeed the new triplet of 
vector bosons to be almost fermiophobic. Hence, only di-boson production, vector boson fusion and triple gauge 
boson production processes can be used to test these models. All these channels require high energy and luminosity 
and will be proper for a future upgrade of the LHC fr3ll - l63j . 

In the strongly coupled scenario, the four site Higgsless model represents a novelty in this respe ct IT3 , l37l I3H] . Its 
phenomenological consequences are quite similar to those of the Minimal Walking Tcchnicolour [35l . |36| . The four site 
model, based on the 

SU(2) L x SU(2) 1 x SU{2) 2 x U(1) Y (111.16) 

gauge symmetry, predicts two neutral and four charged extra gauge bosons, Z\^. and WfHjj and is capable to satisfy 
EWPT constraints without necessarily having fermiophobic resonances. Within this framework, the more promising 
Drcll-Yan processes become particularly relevant for the extra gauge boson search at the LHC. 

The four site Higgslcss model is described by four free parameters: g^ v , g^y, Mzi, Mz2 that is the two vector 
couplings between Zi^-bosons and SM electrons and the two Z12 masses (charged and neutral gauge bosons are 
degenerate). 
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FIG. 6: Parameter space in the plane (fffvi^fez) where g\ v is the ^2-boson vector coupling to SM electrons, and Mz2 is the 
^2-boson mass. The red solid lines restrict the area allowed by EWPT and unitarity. One sample case has been considered: 
2=0.8. 



In terms of the mass eigenstates, the Lagrangian describing the neutral current interaction is given by 3 



(111.17) 



where g x 2V , g( 2A are the vector and axial couplings of the extra Z\i gauge bosons to ordinary matter. In the above 



formula, we have included the g' coupling in the definition of g[ 2V and g[ 2A . 

The energy range, where the perturbative regime is still valid is plotted in the left panel of Figf5]for different values 
of the ratio z = Mzi/Mz2- Owing to the exchange of the extra gauge bosons, the perturbative unitarity violation can 
be delayed up to an energy scale of about y/s ~ 3 TeV. Hence, the mass spectrum of the new particles is constrained 
to be within a few TeV. 

In the past, the only way to combine the need of relatively low mass extra gauge bosons with EWPT was to impose 
the new particles to be fermiophobic. In the four-site Higgsless model, this strong assumption is not necessary anymore. 
In the right panel of Fig|5j we plot the bounds on the vector couplings of the Zi^-bosons to SM electrons coming 
from the EWPT expressed in terms of the ei ; 3 parameters (e2 is uneffective due to SU(2)-custodial symmetry). 
The outcome is that £3 constraints the relation between the two couplings, while e\ limits their magnitude. 

allows to reconcile unitarity and EWPT bounds, leaving a calculable and not fine-tuned parameter space, where 
the new gauge bosons are not fermiophobic. 

Using the linear relation shown in the right panel of Fig|5l we can express the Zi-boson vector coupling to SM 
electrons as a function of the Z 2 -boson vector coupling to SM electrons (by computing back the bare-parameters of 
the Lagrangian, all other Zi^-boson-fermion couplings can be simultaneously derived). In this way, the number of 
independent free parameters describing the four-site Higgsless model gets reduced to three. We choose the following 
physical observables: Mzi, g^vi an d z - I n terms of these new variables, the parameter space allowed by EWPT and 
perturbative unitarity is shown in Fig|6]for one representative z- value: z=0.8. The outcome is that one can reconcile 
unitarity and EWPT bounds, leaving a calculable and not fine-tuned parameter space, where the new gauge bosons 
are not fermiophobic. 

Compared to the popular extra U'(l) theories summarized in TablcU the four-site Higgsless model does not predict 
fixed values for the couplings of the extra gauge bosons to ordinary matter. One has indeed a parameter space 



3 For details see [H H3] . 
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TABLE II: Model predictions and current constraints. The direct limits above on the Z' mass, M z , , are the result of the analysis 
performed in this paper while the best indirect limits, M\, , come from either electroweak (e) fits or contact (c) interactions at 
LEP2 0. 

bounded but still large enough to accomodate rather sizeable Z2-boson couplings to SM fermions. They can range 
from zero to the order of SM couplings. We use this framework in the next sections, when discussing four-site model 
properties and limits. 

IV. APPLICATION OF MODEL INDEPENDENT APPROACH TO THE BENCHMARK MODELS 

A. Current limits from Tevatron 

As discussed above, collider limits on the Z'-boson mass can be presented via contours in the c u — q plane, with 
every contour corresponding to a well defined Mz>- Simultaneously, in the same plane, one can also show the values 
of Cd lU couplings allowed by a specific Z' model. As a result, one can immediately visualize and derive a mass bound 
on the Z'-boson predicted by that particular model. 

We start presenting our results at the Tevatron which is running at yfs = 1.96 TeV. We use the most recent 95% 
C.L. upper bound on a(pp — >• Z' — >• £ + £~) reported at the ICHEP 2010 conference by the DO collaboration for the 
di-electron channel [HI, [66[ where the AM/Mz> — 15% cut was used in the analysis. This limit is shown in FigJT] (left 
panel), together with its 'translation' into the c u — Cd plane for different Mz> masses as shown in Fig J7] (right panel). 

In the following, we use FigJT] (right panel) to interpret current limits from Tevatron, and to derive mass bounds 
on the Z' boson predicted in the classes of models described in the previous section. The results are shown in Fig. [5] 
The top-left panel displays the contour representing E§ Models in the c u — Cd plane, the top-right panel shows the 
Generalised Left-Right Models (GLR), the bottom- left panel contains the Generalised Sequential Models (GSM), 
and finally the bottom-right one gives the Four Site Higgslcss Model (4S). In the first three mentioned panels, the 
colour code corresponds to four equidistant intervals for the mixing angle in the [— 7r/2,7r/2] range for _Eg, GLR, 
GSM models, represented by continuous and closed contours. The black dots on these contours denote the popular 
benchmark models quoted in Table |U In the bottom-right panel, which shows the parameter space of the Four Site 
Model, the colour indicates different mass values for Mz 1 and Mz 2 - The line style distinguishes the Z\ mass (solid 
line) from the Z 2 mass (dashed line). For the Z\ boson, the following mass values have been chosen: Mzi= 480 (red), 
800 (green) and 1600 (blue) GcV. For the chosen sample of free parameters, z = 0.8, the corresponding values for the 
Z-i bosons are: Mzi— 600, 1000 and 2000 GcV shown with the same colour coding. 

Several comments are in order. The first remarkable fact is that there is almost no overlap between contours for 
Eq, GLR, GSM models. This means that, if a Z 1 boson will be discovered and its cross section will be measured with 
a reasonable accuracy, these classes of Z 1 models can be well distinguished using just this basic information. The 
second remark concerns the experimental sensitivity to the Z' production within different models. Comparing the 
four plots, it is clear that the highest sensitivity is to the GSM class of models. In particular, the Q-model can be 
excluded at the Tevatron up to masses Mz 1 > 1260 GeV. Among the GLR models, which are second in terms of the 
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FIG. 7: Left: 95% C.L. upper bound on cr(pp — ► Z' — > e + e _ ) obtained at the Tevatron with integrated luminosity L=5.4fb _1 
by the DO collaboration. Also shown, Z' mass bounds within the extra U(l)' models displayed in the legend. Right: same 95% 
C.L. upper bound on o(pp — > Z' — > e + e~) as above, but translated in the c u — Cd plane into contours corresponding to different 
Mz> masses. The low mass models refer to the Stueckelberg Extensions of the SM (StSM) discussed in [tm l6fi|. 



experimental sensitivity to a Z' boson, the Y-model can be already excluded up to Mz' > 1125 GeV with the current 
Tevatron data. Interestingly, the lowest experimental sensitivity is to E§ models, that is one of the most popular 
class of Z' models. Within this class, the strongest limit can be derived for 8 £ [— tt/4, 0] providing the mass bound 
Mz> > 955, as one can read from the red-coloured part of the E§ contour. 

The 4S class of Z' models must be considered separately. First of all, it predicts two Z 1 bosons with two different 
masses. Secondly, the parameter space of the 4S model is described by more than just one parameter, so the model 
would be represented by an area in the c u — Cd plane rather than by a contour. In order to interpret Fig(8] (bottom- 
right panel) and following analogous figures correctly, a clarification is needed. In the 4S model, the two extra gauge 
bosons can decay into both SM fcrmion pairs and boson pairs. While the contribution to the total width coming from 
the decay into fermion pairs is linear in the extra gauge boson mass, the contribution from the diboson decay grows 
with the third power of the extra gauge boson mass (see Eq. III.6[) . As a consequence, and oppositely to the other 
perturbative gauge models, the Z\^ boson branching ratio into lepton pairs acquires a mass dependence (see Sec. Ill Bl 
for details). This reflects into a mass dependence of the c u and Cd parameters which parametrize the 4S model. So, 
FigEl (bottom-right) should be interpreted as the full parameter space of the four site model projected into the c u — Cd 
plane. To simplify the visualisation of this area, we have varied the vector coupling between the Z2 boson and SM 
electrons, g^y, within the allowed region of Fig. [5] for the sample scenarios: z=0.8 and Mz2= 600, 1000 and 2000 
GeV. This setup should give a full representation of the parameter space, Mz2= 600 GeV being the minimum allowed 
mass and Mzi — 2000 GeV being close to the maximum value of the mass permitted by unitarity. The parameter 
space for these values of M Z 2 an d the respective M Z \ = 0.8 x M Z 2 is presented in Fig. IS^bottom-right) by coloured 
lines (see caption). Whenever the coloured line describing a given Mzi,zi value for the 4S model crosses the black 
contour corresponding to the same mass value, that crossing point would give the experimental sensitivity to a Z\^ 
boson with that mass. The portion of the coloured line above this crossing point would represent the excluded region 
in the c„ — Cd parameter space. To clarify the interpretation, let us consider the following examples. For Mz 2 = 1000 
GeV and Mz 1 = 800 GeV (green lines), one can see that the parameter space for the Z\ boson (solid green line) is 
greatly excluded in the region of the c u — Cd plane above the black contour line labelled by 800 GeV. Only a small 
portion of the c u — Cd parameter space is instead excluded for the Z2 boson (dashed green line). The c Ut d couplings 
are indeed bounded to be c Ui d < 10~ 3 , as one can see from the crossing point between the dashed green line and the 
black countour labelled by 1000 GeV. Since, the two extra gauge bosons would be simultaneously produced, from 
the discussed green lines one should deduce that the most restrictive bound on the c U) d couplings comes from the Z\ 
boson. If not observed, the crossing point of the solid green line representing Mzi = 800 GeV with the black contour 
at 800 GeV would give the bound: c Ut d < 210 -4 . Consider now a second scenario: Mzi = 2000 GeV and Mzi — 1600 
GeV. This is represented by blue lines. In this case, the Tevatron has no sensitivity at all to the two extra gauge 
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FIG. 8: 95%CL limits on M z > in the c u — Cd plane based on the 2010 analysis of the di-electron channel performed by the DO 
collaboration at L=5.4fb~ 1 . Top-left, top-right, bottom-left, bottom-right panels present results for Ea Models, Generalised 
Left-Right Models (GLR), Generalised Sequential Models (GSM) and Four Site Higgsless Model (4S), respectively. The colour 
code corresponds to four equidistant intervals for the mixing angle within the [—n/2,n/2] range for Ea, GLR, GSM models. 
The black dots on these contours correspond to the benchmark models listed in Tabled] F° r the 4S model, the colour denotes 
different values for Mz x and Mz 2 - The line style distinguishes the Z\ mass (solid line) from the Z2 mass (dashed line). For 
the Z\ boson, the following mass values have been chosen: Mzi~ 480 (red), 800 (green) and 1600 (blue) GeV. For the chosen 
sample of free parameters, z — 0.8, the corresponding values for the Z2 bosons are: Mzi= 600, 1000 and 2000 GeV shown with 
the same colour coding. 



bosons. The solid blue line representing Mzi — 1600 GeV, in fact, never crosses the corresponding black contour line 
labelled by 1600 GeV. And the same is true for the dashed blue line representing Mzi = 2000 GeV. 

Even in more complicated multi-resonance scenarios, as in the 4S model, the c u ^d representation allows one to 
visualize directly via the black contour lines, up to what mass value the experiment could be sensitive. If no signal of 
new physics is observed, the bound on the mass can be translated into limits on the c u _d coefficients. From there, one 
can then trace back exactly the lagrangian parameters of the 4S model which are in turn excluded. In case of new 
physics discovery, the c u ^d approach allows one to uniquely determine the c u _d values corresponding to that observed 
mass. Since the c Ut d coefficients arc strictly related to the new gauge boson couplings, this in turn enables one to 
extract informations on their size. In the 4S model, the coupling between the Zi boson and SM electrons, g^y-, grows 
linearly with c u and Cd as shown in Fig. [§1 Any mass measurement therefore translates into a coupling determination. 
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Moreover, g^y is one of the three free parameters of the model. The measurement of the mass of the new gauge 
bosons would therefore allow one to derive direct informations on the bare lagrangian parameters. 




FIG. 9: g2v fermion-boson coupling versus c tl (left) or Cd (right) in four-site model 




FIG. 10: Left: CMS limit on the Z' boson production cross section in the di-electron channel, normalized to the SM Z boson 
cross section, as a function of the Z' mass. The limit is projected at 500pb _1 [glj. Right: Limits in the c u — Cd plane, based 
on the projected LHC 500pb _1 limit shown in the right panel. In the c u — Cd representation, the limits appear as contour lines 
corresponding to different M z > values. 



B. Expected LHC potential at 7 TeV to probe Z' models 

We now explore the LHC@7TeV potential to test the classes of Z' models under consideration. We use the projected 
limits from LHC. In particular, we rely on the limits given by the CMS Exotica group for 500pb _1 of integrated 
luminosity which hopefully will be available in about one year from now. This limit is shown in Fig llOl (left panel), 
which is taken from the public web page of the CMS Exotica group (69|. The projected 500pb _1 limit from CMS is 
given as a ratio oz' /&z, where oz is the Z-boson production cross section in the 60 < m ee < 120 GeV window and 
we have converted this limit into the limit on the NNLO production cross section for the Z' boson shown in terms 
of c Uf d coefficients in c u — Cd plane for different Z' masses given in the right panel of Fig llOl This representation is 
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the analog of what done before at the Tevatron. Comparing FigJT] and FigflUl one can observe the strong gain in 
sensitivity one gets at the LHC@7TeV with 500pb -1 , compared to the Tevatron with 5.4fb -1 , at fixed c Ut d value. 





FIG. 11: Limits in c u — Cd plane are based on projected CMS 500pb 1 limit. 95%CL limits on M Z ' in c u — Cd plane based on 
on projected CMS 500pb _1 analysis of di-electron channel. The legend scheme is the same as in Fig. [H] 

Now we can estimate the LHC@7TeV potential for deriving bounds on Z' models at 500pb _1 . The results are 
shown in Fig. [TTJ where the legend scheme is the same as in Fig. [5] From Fig. [TT] and Fig. [HJ one can see that for the 
models with small-intermediate values of Z' boson couplings to SM fcrmions (that is Eq models, GLR-modcls partly, 
and some GSM models), the LHC@7TcV can extend the limits on Mz> by about 500 GeV when compared to the 
Tevatron. For example, the limit on the SM-likc Z' boson could be extended from 1020 GeV to about 1520 GeV. 
On the other hand, the limits for larger c u ,d coefficients and respectively larger masses could be extended in the near 
future up to a 2 TeV scale, which is unreachable at the Tevatron. For the Q-model, belonging to the GSM class, the 
mass bound could be improved from 1210 GeV (current Tevatron) to 2250 GeV. Regarding the 4S model, one can 
see that the scenario characterized by Mzi = 800 GeV and Mzi = 1000 GeV could be totally excluded in the c u — Cd 
plane shown in Fig. 1111 The solid green line, representing the Z\ boson parameter space, lies in fact beyond the 800 
GeV black contour line in the dispayed plane. No improvement, compared to the current Tevatron, would instead be 
possible for the scenario: M Zl = 1600 GeV and M z . 2 = 2000. The LHC@7TeV and 500 pb" 1 will not have sensitivity 
enough, as one can deduce from observing that the blue lines never cross the black contour lines labelled by their 
respective mass values. Exploring this region of the 4S parameter space would require higher integrated luminosity 



19 



and preferably higher collider energy. 





FIG. 12: 5a Discovery limits in the c u — Cd plane, based on the 500pb -1 projected analysis of the Z' production in the di- 
electron channel performed by the CMS collaboration. The legend scheme is the same as in Fig. [8] The upper(dashed) and 
lower(dot-dashed) contours correspond to the uncertainty in Fig. HOf left) reflected in the width of the band. 

We also use Fig. [10] to estimate the LHC@7TeV discovery limits for 500pb -1 . In this analysis, we assume that 
the significance grows as y/L, where L is the total integrated luminosity and that the signal over background ratio is 
constant for the same selection cuts. The last assumption is motivated by the fact that for a chosen invariant di-lepton 
mass window the qq parton densities are very similar for signal and background processes and defined mainly by the 
y/I value. 

The LHC discovery potential for various Z' models is shown in Fig. [I2]in the c u — Cd plane. The legend scheme is the 
same as in Fig. [5J The upper(dashed) and lower(dot-dashed) contours correspond to the uncertainty in Fig. HoTlcft) 
reflected in the width of the band. One can see that discovery limits are typically 150-200 GeV lower than exclusion 
ones. 
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V. IMPACT OF Z' WIDTH ON SEARCH STRATEGIES 

Invariant mass distributions may be examined in a number of ways for evidence of resonant structures. The 
sensitivity of any particular approach has a dependence on the intrinsic width of any possible resonance. The simplest 
approach is to bin the invariant mass distribution and determine the compatability of the number of events in any 
bin with the Standard Model prediction. A p-value may be used to quantify this compatability. In this approach 
the width of the bins for optimal sensitivity depends on the intrinsic width of possible resonances and the detector 
resolution. In the case where the width of the resonance is much smaller than the detector resolution this parameter 
defines the optimal bin size. For intrinsic widths comprable to the detector resolution then the optimal width depends 
on both of these parameters. 

An alternative is to use a parameterization of the expected distribution in the two alternative hypotheses of, a 
distribution resulting only from Standard Model Physics and one resulting from the addition of a new physics process. 
A comparision of some measure of the quality of the fit in both cases allows a determinatin of the probability of the 
presence of New Physics. In this case the functional form of the resonance structure is typically taken to be some 
variant on a convolution of a Gaussian and a Breit-Wigner. For a resonance where the width of the Breit-Wigner 
is small in comparision to the width of the Gaussian the sensitivity of the search is insensitive to the width of the 
Breit-Wigner. Such circumstances result in the greatest possible experimental sensitivity. For resonances with large 
widths compared to the experimental resolution then the Breit-Wigner width must be included as a further parameter 
in any fits. In cases where the interference effect is large this must be included in the fuctional form used to fit to the 
invariant mass distribution. We show here that in the mass regions which will provide the greatest sensitivity for low 
integrated luminosities at the LHC the interference effect is negligible and may be ignored without compromising the 
search sensitivity. 

The above descriptions of possible methods of searching for a resonance in an invariant mass distribution illustrate 
that a knowledge of the width of the resonances being searched for has an impact on the search procedure used. In 
order to obtain the best sensitivity, it is thus important to have a knowledge of the magnitude of the widths from 
New Physics models. The results of such searches depend on the assumptions made in the analysis and can't be easily 
interpreted in other circumstances. 

Besides that, in the case of a discovery the Z' boson mass and decay width will be determined from fits to the 
reconstructed invariant mass distribution of di-lepton candidate events. In this section, we thus focus on the prediction 
and the possible measurement of the total decay width, comparing the various classes of models under consideration. 

All extra U'(l) theories summarized in Table U make the assumption, for sake of simplicity, that the Z'-boson 
decays purely into SM fermions. Under this approximation, the total decay width is given by Eq |II.7l and its value 
never exceeds a few percent of the corresponding mass (Fz>/Mz> < 3%), as shown in Table [TT1 This property has a 
direct implication on the possibility to measure the Z' decay width at the LHC, being correlated to the experimental 
di-lepton mass resolution. If indeed Tz 1 > R, being R the mass resolution, one can have direct access to the decay 
width of the observed spin-1 particle. During the early stage of the 7 TeV LHC, the expected di-electron mass 
resolution is about Rlhc = 2%M [70l [7lj. As a consequence, within the majority of models summarized in Table 
HT1 the total Z'-boson width is hardly measurable. This is visualized in Fig. [13J where the ratio between Z' width 
and mass is plotted as a function of the mixing angle parametrizing the three classes of models listed in Table Q] Eq, 
GLR and GSM. All Eq inspired models predict a quite narrow Z' boson. For some benchmark model within the 
generalized Left- Right class (GLR), the ratio becomes slightly bigger than the early LHC resolution. The scenario 
changes when we consider generalized sequential models (GSM). Here, the width over mass ratio is well above the 
early 2% resolution of the LHC. Another example of measurable decay width is given by the four-site Higgsless model. 
Here, in fact, in most part of the parameter space Fz 1-2 /Mz 1 2 > 2%. This property is shown in Figfn]ior different 
values of the free z parameter. In this case, the distinctive behaviour is due to the fact that the Zi ;2 -bosons predicted 
by the four-site model decay preferebly into the diboson channel: Z\ — > WW and Zi — > W\W . Thus, their width 
grows with the third power of the corresponding mass, as shown in Eqs lII.12III.13l making it wider compared to the 
other models. This feature is common to all Higgsless and Technicolor models. 

This discussion is appropriate when considering the Z' boson production in the di-electron DY channel. For di-muon 
final states, the analysis would change drastically. The estimated di-muon mass resolution during the early stage of 
the 7 TeV LHC is in fact around R ~ 8% |7l|. Hence, only in a very restricted range of the GSM mixing angle 
the Z' width could be measurable. An exception is given by the four-site Higgsless model, where the Z^ boson width 
could be determined in a large portion of the parameter space. 

The discussed differences between classes of models should be taken into account for improving search strategies 
and possibly measurements. 
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FIG. 13: Z 1 boson width over its mass as a function of the mixing angle parametrizing the GLR and GSM class of models 
given in Table [I] The colour code corresponds to four equidistant intervals in the [—n/2,n/2] region. The black dots on the 
contours correspond to the benchmark models listed in Table [I] 




FIG. 14: Left: Zi-boson width over its mass as a function of Mzi for different values of the free z-parameter. Right: same for 
the Z2-boson. We assume the maximal ^2-boson coupling to SM fermions. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have discussed the prospects for setting limits on or discovering spin-1 Z' bosons using early 
LHC data at 7 TcV. In order to facilitate the connection between experimental data and theoretical models, we have 
advocated the narrow width approximation, in which the leptonic Drell-Yan Z' boson production cross-section only 
depends on the Z' boson mass together with two parameters c u and c^. These variables provide a convenient way of 
expressing the experimental limits or discovery information about the mass and cross-section which enables a direct 
comparison to be made with the predictions arising from theoretical models. The experimental limits on the Z' boson 
cross-section may be expressed as contours in the c u — Cd plane, with a unique contour for each value of Z' boson 
mass. If a discovery is made then the measurement of the mass and cross-section corresponds to some unique contour, 
or in practice a unique band including the error. Such contours may be compared to the theoretical prediction of 
c u and Cd arising from particular models, enabling limits to be set on models or a discovery to discriminate between 
particular models. However the application of this strategy requires the experimental cross-sections to be properly 
defined and the theoretical cross-sections to be accurately calculated as we now discuss. 

On the experimental side, we have seen that the use of the narrow width approximation requires an appropriate 
di-lepton invariant mass window cut around the mass of the Z' boson. The effect of the cuts is rather subtle since 
it depends on both the width of the Z' boson and the energy of the collider, with higher widths and lower collider 
energies leading to more prominent signal tails at low invariant masses. Fortunately we have seen that at LHC energies 
the suitable experimental cut is rather insensitive, especially for models with lower Z' boson widths, and furthermore 
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may be optimised at a unique value suitable for all models, although there is some unavoidable dependence on 
the Z 1 mass. One important conclusion is that whatever cut is chosen should involve invariant masses well above 
100 GcV otherwise interference effects will be dominant. In summary, we have demonstrated that this cut plays a 
crucial role: it diminishes a possibly huge model-dependent interference effect, removes model-dependent shape of the 
M^+i- distributions in the region of low Mg+g- especially for the case of large Z' width and brings into agreement 
NWA and FWA cross sections. On the theoretical side, we have evaluated cross-sections at NNLO using updated 
ZWPROD package. One should stress that K^nlo factors are depend on the Z' mass and we have tabulated them 
for convenience for both the Tevatron and LHC. Moreover the Knnlo ar e very PDF dependent and one should 
specify which PDF is being used to apply a respective K^nlo- 

We have applied the approach above to two quite different types of Z' models: perturbative gauge models and 
strongly coupled models. Among the perturbative gauge models we have studied three classes: E§ models, left-right 
symmetric models and sequential standard models. Each class of model is defined in terms of a continuous mixing 
angle variable. This enabled infinite classes of benchmark models to be defined, rather than just a finite number of 
models, where for each class of model, the respective angles serve to parametrize different orbits in the c u -Cd plane. 
These orbits turn to be non-overlapping for these three classes of models. A limitation of this approach is that it 
ignores the effect of the SUSY and exotics (and right-handed neutrinos) on the width Yz> ■ Assuming only SM particles 
in the final state we have calculated the widths of the benchmark classes of models and seen that the perturbative 
models generally involve relatively narrow widths (which however can be increased if SUSY and exotics are included 
in the decays), while the strongly coupled models involve multiple Z' bosons with rather broad widths. We have 
also commented on the significance of the width on search strategies which if measured would a complementary for 
a discrimination of the underlying Z' model as well as it would allow to test non-SM Z' decays menationed above. 
Another limitation of our approach is that it ignores the effects of Z — Z' mixing which is quite model dependent. 
However such effects must be small due to the constraints from electroweak precision measurements, so such effects 
will not have a major effect on direct collider searches considered here, although of course they will affect the precise 
vector and axial couplings (see for example [(111] where the U(1)n vector and axial vector couplings are calculated 
including the mixing effects). As regards the strongly coupled models, we only considered the four-site Higgless model 
which contains just two excitations of the Z (and W) bosons. However it is representative of models of walking 
technicolour models and the KK excitations of the Z(and W) which is considered for the first time in the c u -Cd plane. 
It is clearly seen that the associated Z' bosons may easily be distinguished from those of the perturbative gauge 
models. 

In conclusion, our results support the use of the narrow width approximation in which the leptonic Drell-Yan Z' 
boson production cross-section only depends on the Z' boson mass together with two parameters c u and Cd- However, 
as discussed in this paper, care must be taken concerning the experimental cuts and the theoretical Knnlo factors 
tabulated here must be included correctly. Providing the experimental cross-section is appropriately defined, according 
to the recipe we provide in Fig|3J and the theoretical cross-sections are properly calculated at NNLO, we have shown 
that such a strategy is safe, convenient and provides the most unbiased way of comparing experiment to theoretical 
models which avoids any built-in model dependent assumptions. The experimental limits or discovery bands may then 
be reliably confronted with the theoretical predictions on the c u -Cd plane as shown in Fig |8llllI2[ which represent 
the main results of our study, leading to the new limits which we derive here for the Tevatron and to the projected 
limits for LHC. The results show that the LHC at 7 TeV with as little data as 500 pb _1 can either greatly improve on 
current Tevatron mass limits, or discover a Z' , with a measurement of the mass and cross-section providing a powerful 
discriminator between the benchmark models using this approach. 
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Appendix A: ifjviVLO-factors and cross sections for the pp(p) — > Z' + X production process at the Tevatron and 

the LHC@7TeV 
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TABLE III: Values of TfjvjvLo-factors and cross sections for the pp(p) — > Z' + X production process at the Tevatron and the 
LHC@7TeV, obtained with ZWPROD package for MSTW08 PDF. 
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TABLE IV: Values of TfjviVLO-factors and cross sections for the pp(p) — > Z' + X production process at 
LHC@7TeV, obtained with ZWPROD package for CTEQ6.6 PDF. 
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